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Abstract 

We present an algorithm to generate positive braids of a given length 
as words in Artin generators with a uniform probability. The complexity 
of this algorithm is polynomial in the number of strands and in the length 
of the generated braids. 

As a byproduct, we describe a finite state automaton accepting the 
language of lexicographically minimal representatives of positive braids 
that has the minimal possible number of states, and we prove that its 
number of states is exponential in the number of strands. 



1 Introduction 

Recently, many research papers have appeared that describe effective compu- 
tations with braids on n strands. In many cases, the results of performing a 
certain computation for a set of 'random' or 'pseudo-random' braids of a given 
length are reported. Usually, the authors generate so-called positive braids, that 
is, elements of the submonoid of the braid group B n generated (as a monoid) 
by the standard Artin generators {ai, . . . ,<7„_i}. Methods commonly used to 
generate such positive braids are the following: 

1. In order to generate a positive braid of length fc, choose k times an element 
of {ai, . . . , <7 n _i} with uniform probability, and form their product. 

2. In order to generate a positive braid of canonical length k, choose k times 
a simple braid with uniform probability on the (finite) set of simple braids 
and compute their product; if the obtained braid has canonical length 
smaller than k, discard it and try again. 

However, none of these procedures generates braids with a uniform distribution; 
some braids are far more likely to appear than others. For instance, if one uses 
the first procedure to produce positive braids of length 6 in B4, the probability 
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of obtaining the braid (cti) 6 is 3 -6 , as there is only one way to write this braid 
as a positive word in Artin generators. On the other hand, the probability of 
obtaining the braid A = o\020~\o~-so~io~\ is 16 ■ 3 -6 , as there are 16 distinct ways 
to write A S B4 as a positive word in Artin generators [5]. That is, the above 
procedure is 16 times more likely to generate A than to generate (cti) 6 - This 
bias becomes more dramatic as n and the length of the braids involved increase. 

In this paper we shall give a procedure to generate random positive braids of 
given length with a uniform probability. That is, given n and k as input, the 
algorithm generates a positive braid in B n of length k, in such a way that the 
probability of obtaining any given braid is l/x ni k, where x n ^ is the number of 
positive braids of length k in B n . 

The structure of the paper is as follows. Section [2] describes the basic idea 
of generating uniformly random positive braids via lexicographically minimal 
representative words (lex-representatives, for short). In Section [3l we count the 
braids in of a given length k using a result by Bronfman. In Section |4l 
we develop a description of lex- representatives, which will be used in Section [S] 
to count the lex-representatives that start with a given prefix, completing the 
description of our algorithm. 

It is known that the set of lex- representatives is a regular language. In Section^ 
we show that our description of lex-representatives yields an acceptor for this 
regular language that has the minimal number of states, and that the number 
of states is exponential in n. This shows, in particular, that using standard 
language theoretical techniques to generate uniformly random positive braids is 
not efficient. Finally, in Section we analyse the complexity of our algorithm 
and give timing results. 

Acknowledgements: We thank Pascal Weil, Frederique Bassino and Cyril 
Nicaud, for suggesting to us the use of automatic structures to generate random 
elements, which led us to realise that our sets of minimal forbidden prefixes 
provide a minimal finite state automaton. 

2 Structure of the algorithm 

The basic idea of our algorithm to produce random positive braids is the fol- 
lowing. Let An = {(7i, . . . , cr„-i} be the set of standard generators (or atoms) 
in £?+. Let A n be the free monoid generated by A n . We know there is a mor- 
phism of monoids b : A* n — > B+ which sends each element in A* n (a word in 
<7i, . . . , <7 n _i) to the positive braid it represents. As we saw above, the map b 
is not injective, so we are going to define a section of b. For that purpose, we 
will order the elements in A* n having the same length using <i ox , which is the 
lexicographical order in which o\ < &i < ■ ■ ■ < <x„_i. 

Definition 2.1. Given (3 £ we define the lex-representative of (3 to be 

u>(fi) =min{b-\f3)} EA* n . 
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In other words, w(/3) is the smallest positive word, with respect to <i ox , that 
represents (3. Notice that this is well defined as is a homogeneous monoid 
(all words representing a given element have the same length). For instance 
w(<73<7i) = cti<73, uj{<J2<Ji<y2) = o\02(J\ and for A G B4, w(A) = o\02<J\o-a<J2<J\- 
Here the arguments of to are braids in _B+ , whereas the images are words in A* n . 
It is clear that lu : B+ — > A* n is a section of b, so it is injective. Moreover, 
as -B+ is homogeneous, we can define {B+)j~ to be the set of positive braids of 
length k, and L n: k to be the set of lex-representatives of length k (L stands for 
language). Then lu : (B^)k — * ^n,fc is a bijection, whose inverse is b\ L . 

Our algorithm will produce a random element of L n ^ with uniform probability. 
By the above arguments, this is equivalent to producing a random positive braid 
of length k with uniform probability on (B^)k- There are three main steps: 

1. Determine the size of L n _k, that is, compute x n< k = | (£?+)& |. 

2. Choose a random integer r between 1 and x„.fc. 

3. Find the r-th word in L n ^, where L n ^ is ordered by <i cx - 

The second step poses no difficulty (assuming we know how to generate random 
integers), so our task is to be able to perform steps 1 and 3 in polynomial time. 
It is useful to identify the lex-representatives of braids in B+ as the vertices 
of a rooted tree, with the elements of L n ^ being the vertices at depth k and 
edges given by the prefix partial order in the monoid A* n . (That is, the root 
corresponds to the trivial word e, and given a word w G L n , s and <ii G A n such 
that woi € L n s+ i, there is an edge labelled Oi joining w to woi ) 
Figure [1] shows the tree of lex- representatives in B\ truncated at depth 3; 
the elements of L4.3 correspond to the leaves of the truncated tree. We have 
I Z/4,0 1 = lj 1-^4. i| = 3, I Z/4 5 2 1 = 8 and ^4,3! = 19. Notice that there are no 
vertices corresponding to the words a^<J\ and cr2CiC2, as these words are not 
lex-representatives. The fact that there are 19 leaves means that there are 
19 positive braids in B^ of length 3. Hence, in this example, our algorithm 
will choose a random number r between 1 and 19, and it will look for the 
r-th leaf of the rooted tree, counting from left to right, as leaves are ordered 
lexicographically by construction. 

Of course, it is not efficient at all to try to compute the whole set L n ^k, as it has 
exponential size with respect to k. We shall overcome this problem by counting 
suitable subsets of leaves of the truncated rooted tree, as follows. 

Definition 2.2. Given a word w G A* n and an integer m £ {0, . . . ,n — 1}, we 

define x n ^k{w,m) to be the number of words in L n ^ of the form ww' , where 
w' G does not start with o\ , . . . , a m . 

Obviously, x n ^(w,m) — if w is not a lex-representative, or if \w\ > k. More- 
over, if w is a lex-representative, then x rit k{w,n — 1) = 1 if k = \w\ and 
x n ,k{w,n — 1) = otherwise. We will show in Section [3] how to compute 
£n,fc(e, 0) = Xn t k, and in Section [6] how to compute a; ni fe(e,m) for m > 1, and 
also x nt k{w,m) where w = w'o~j and m > j — 1, in polynomial time and space 
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Figure 1: The rooted tree corresponding to i 4j3 



with respect to n and k. If we assume this to be known, we can explain our 
main algorithm with some more detail: 

Theorem 2.3. Suppose the time and space required to compute x ni k(w,m), 
where m > if w — e and m > j — 1 if w = w' ' o~j , is polynomial in n and k. 
Then there is an algorithm generating a random element of with uniform 

probability, whose time and space complexity is polynomial in n and k. 

Proof. In order to generate a random element of (_B^) fe , we compute x n .k and 
choose a random integer r between 1 and x n .k- It just remains to determine the 
r-th element u>( r ) of the chain L n ^ with respect to <i cx . We will find the word 
W( r ) letter by letter. (Pseudocode is given in Algorithm [T]) 

Suppose that we have computed a prefix w of u>( r ) and know that there are v 
words in L n ^ of the form to', where w' € A* n and u>( r ) <iox ww' . (Initially, 
w = e and v = x n .k — r.) If w ^ e, let <jj be the last letter of w; otherwise 
let j = 0. Let o~i be the next letter of W( r y Observe that i > j — 1, since 
otherwise W( r ) would not be lexicographically minimal. By Definition 12. 2\ we 
have to > i if and only if x n ^(w,m) < v. Hence, we have i = minS 1 , where 
S = {to G {j — 1, . . . , n — 1} I x n _k{w, m) < ^} 3 n — 1. We can determine 
i using binary search in at most [log 2 (n — 1)] steps, each step requiring the 
computation of one number of the form x n ^{w, to). We then replace w by wa% 
and v by v — x 7h k(w, «), and proceed to the next letter of W( r ). (While this does 
not affect the complexity, we remark that x n ^{w,i) is known from the binary 
search, so updating v just requires one subtraction.) 

As u>( r ) has length fc, at most fc|~log 2 (n — 1)] steps are required to find W( r y, 
since each step is polynomial in n and fc, the claim follows. □ 

Example 2.4. Suppose we want to generate a random braid of length 3 in B\ 
(see Figure^). Compute 1:4,3 = 19. Choose a random number r between 1 and 
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19. Say r — 16, so v — 3. Now compute 2:4,3(6, 2) = 4 > 3 = v . Thus, the first 
letter ofwn^ is 03. Since 2:4,3(6, 3) = 0, the value of is remains unchanged. 
Now we find the second letter ofwn&s, which must belong to the set {(72, 03} (as 
the first letter is 03/ We compute 2:4,3(03,2) = 2 < 3 = v. Hence, the second 
letter of Wn^ is 02, and we must replace v by v — 2:4,3 (03, 2) = 1. 
Finally, as the second letter o/w(i6) is o~2, the third one could be any letter in 
{<7i, 02, 03}. We compute a;4 1 3(o'3<72,2) = < 1 = v, then 2:4,3 (0302 } 1) = 1 < 
1 = v, and conclude that the third letter of iU(i6) is o~\. Thus, whq) — o^o^ci- 

Algorithm [1] gives pseudocode for finding wt r \, assuming we know how to com- 
pute x n> k{w, m) as in the statement of Theorem 12.31 It contains references to 
two algorithms which will be introduced in Section [3] respectively Section [5] In 
Section [3] we will use a known result [2] to compute x„ i fc(e,0) = 2; n! fc (line Q] 
of Algorithm [T]) . An algorithm for computing all other required instances of 
Xn,k{w,m) (line [8] of Algorithm [lj will be given in Section [6] 



Algorithm 1 Producing a uniformly random braid in of length k 
Input: Integers n > 2 and k > 0. 
Output: A braid in B„ of length k. 

1: Compute x n ,k- [Section [3] equations ([2]) and <[3j ] 

2: Choose a random integer r £ {1, . . . , x n ,k} with uniform probability. 

3: w := e ; v := ~ r ; a := 

4: for i := 1 to do 

5: a := max{a — 1, 1} ; 6 := n — 1 ; yit := /* fj, — x n ,k(w, b) at all times */ 
6: while a < b do 
7: m := L^J 

8: Compute x n ,k(w, m). [Algorithm [2] 

9: if x nt k(w,m) < v then 
10: b:=m; fi := x n ^(w,m) 

11: else 
12: a := m + 1 

13: u> := wa a ; f := v — \i /* at this point a — b, hence /i = x n ,k(w, a) */ 
14: return b(w) 



3 Counting all positive braids of a given length 

In this section we use a formula given by Bronfman [5] to describe how to 
compute the number of elements of [B^j, in time respectively space that is 
polynomial in n and k. To our knowledge, the paper [2] was never published, 
but one can find proofs of its main result in p] and also in [7] . Bronfman gave a 
recurrence relation for the growth function of the monoid B+. More precisely, if 
x n ,k is the number of elements of length k in B£, then the formal power series 

G n (t) — S ' x Ut k t k 

k>0 
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is the growth function of B+. Deligne [I] showed that this function was rational, 
namely the inverse of a polynomial H n (t). The following recurrence relation to 
compute this polynomial was given in [2J: 

n 

H n (t) = J2(- 1 Y +1 t® H n-i(t), (1) 
z=l 

where Ho(t) = H\(t) = 1. In particular, H n (t) is a polynomial of degree (™) 
such that H n (0) = 1. Denoting the coefficient of t m in ffi(i) by hi >m , it follows 
from G n (t)H n (t) = 1 that for any j > one has 

2n,j = -(Xn,j-l/ln,l + ^n,j-2^n,2 H h X nfl h n ^) , (2) 

where the sum has at most min{j, f ™J } terms, since h„ im = if m > (™) • 
Therefore, we can iteratively compute the numbers x n $, x ni i, . . . , x n> k, knowing 
the coefficients hni, . . . , /i n ,fe- From (|T|), the latter can in turn be calculated 
with the recurrence relation 

m 

*W = E(- 1 ) <+1 '»m-i I i-G)> ( 3 ) 
i=l 

where we define h m _ i:j _n\ = if j — (2) < 0. Note that in order to compute 
x n> k, only the coefficients h m j for < m < n and < j < h are required. In 
particular, time and space required to compute x n ^ are polynomial in k and n. 

The coefficients h m _j as well as the integers x m j for < m < n and < j < k 
will be used later in our algorithm generating random elements of Note 
that coefficients h m j that have already been computed do not change if n and 
k are changed; one merely might have to compute additional coefficients. Thus, 
any values of h m j and x m j that are computed can be stored once and for all; 
their computation should be thought of as a precomputation. 

4 Forbidden prefixes 

It only remains to compute x ni k(w, m), where m > 1 if w = e, and m > i — 1 
if w = w'oi. (Recall Definition 12.21 ) This calculation, which will be done in 
Section relies on a detailed description of the words w' that can be appended 
to w so that the concatenation ww' is still a lex-representative in terms of what 
we call forbidden prefixes. 
Denote by 

L n = |J L n<k = |J w((B+) fc ) C A* n 

k>0 k>0 

the set of words in {ox, . . . , <r n _i} that are lex-representatives. To simplify 
notation, let M = and M]~ = (B£)k for k £ N. We denote the length of a 
word w £ A* n by \w\, and the length of a braid x £ M by 
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Definition 4.1. For w £ L n , we define the set of forbidden prefixes after w 

as F n (w) = {a £ M | wui(a) ^ L n }. 

For a braid a e M, we define the set of multiples of a as aM = {a/3 | /? S M}. 
The following is an important property of the sets of forbidden prefixes. 

Lemma 4.2. Let w £ L n . If a £ F n (w), then aM C F n (w). 

Proof. For u £ JC n and fceN, let it| fc denote the initial subword of length k of u. 
As a £ F n (w), we have u>(b(w)a) <i cx u^ta), and thus w(&(iu)a)| |iu| <i ex 
For any f3 £ M, the latter implies uj(b(w)a/3) <i ox Lu(b(w)a)uj(P) <i cx ww(a/3), 
proving the claim. □ 

As M is a cancellative monoid, it admits a well-defined partial order where 
a =4 b if there exists c S M such that ac = b\ in this case we say that a is a 
prefix of £>. It is well known that (M, =4) is Noetherian, that is, that there are no 
infinite descending chains with respect to =<; in M. We can then define a special 
subset of the forbidden prefixes. 

Definition 4.3. Form £ L n , the set of minimal forbidden prefixes after w, 

denoted F™ ln (w), is the set of minimal elements, with respect to ^4, in F n (w). 

From Lemma and the Noetherianity of M, one has F n (w) = {J a< =F^( w ) otM. 

The rest of this section will be devoted to the study of F™ ln (w). We will show 
that it is a finite set, we will describe its elements, and we will show how to 
compute them in linear time with respect to n and \w\. 

First we recall a particularly useful property of the prefix order =^ of M : It is 
a lattice order [5], that is, for any a,b £ M there are a greatest common prefix 
a A b and a least common multiple aVt with respect to . In the particular 
case of Artin generators, we have o~i V <jj = oia^ = o~jO~i if \i — j\ > 1, and 
cri V o~j — OiOjOi = o~jO~iOj if \i — j\ = 1. 

Given two elements a,b £ M, there are positive elements a\b and b\a, such that 
a V b = a(a\b) = b(b\a). As the monoid M embeds in the group B n , we can 
write a\b = a _1 (aV6) and b\a = b~ 1 (a\/b). Notice that o~j\ai — cr,; if \i — j\ > 1 
and o-j\ai = a^Uj if \i — j\ = 1. 
We shall frequently use the following property: 

Lemma 4.4. // a, 6, c £ M , then a ^ be is equivalent to b\a =^ c. 

Proof. As b =4 be, it follows that a =^ be if and only if a V b =^ be; the latter is 
equivalent to b^ 1 (a V b) =<; c, so the claim is shown. □ 

The following two results form the basis for an inductive description of F™ ln (ui). 
Proposition 4.5. 

(a) F n (e) = F;r(e) = <t). 



7 



(b) F 7 f»(<j J ) = {a 1 ,...,a^ 2 ,a^ 1 a 3 }. 

(Here and in the sequel, we use the convention that words involving indices 
less than 1 are ignored. That is, F™ la (o~i) — and F£ ral (o~2) — {o~io~2}-) 

Proof. It is obvious that {a%, . . . , <7j-2, o~j-i°~j} ^= Fn{°~j)- Indeed, u\, . . . , Oj-2 
are atoms and ^ F n (o-j), so {a%, . . . , <7j_2, ^ F™ ln (°~j)- 

Conversely, if a £ F n (o~j) then ui(o~jCt) <i ox o-jUj{a), and hence uj^aja)^ <i cx o~j. 
This means that cr^ =^ o^a. for some i < j, whence a admits a prefix from the 
set {a 3 \ai | i = 1, . . . , j - 1} = {o"i, . . . , <7,-_2, c^-icr,-}. 

This proves Claim (0). Claim (Jaj) is trivial. □ 

Proposition 4.6. Let w G L„ wzi/i |w| > 1. If w = vaj, then F^ lln (w) is the 
set of =4-minimal elements in 

{a 1 ,...,a j _ 2 ,a^ 1( r j }U{a j \l3 | /3 G F„ mi »}. 

Proof. The case |w| = 1 is clear from Proposition 14. 51 as F™ m (e) = 0. Suppose 
that \w\ > 1 and let a be a forbidden prefix after w — vaj, that is, vo~jtx){a) L n . 
We will distinguish two cases: either <jjUi(a) is a lex-representative, or not. 
If o-jLo(a) is not a lex-representative, then a is a forbidden prefix after o~j. By 
Proposition 14.51 in this case a admits a prefix a' £ {cri, . . . , o~j^2, a j-i (7 j}- 
If <7jU>(a) is a lex-representative, then o^a is a forbidden prefix after v, hence 
there is (3 G F™ m (w) such that /3 =<; o^a, that is, erj\/3 a. 
We have then shown that if a is a forbidden prefix after w, then a admits a prefix 
at £ {ax, . . . , aj-2,o~j-iC~j} U {o"j\/3 | /3 £ F™ m (w)}. Since the elements of the 
latter set arc forbidden prefixes after w by construction, the claim follows. □ 

The following lemma will give us control over the elements of the form o~j\/3. 
Lemma 4.7. For i G {2, . . . , n — 1} and j G {1, . . . , n — 1}, one has: 



(cri_icr i )(cri_20'i_i) if j = i - 2 

<7i if j = i-l 

o~i-io~iO-i+\ ifj = i + l 

^ai-iO~i otherwise 



For j £ {1, . . . ,n — 1} and l<m<z<n— 1, one /las." 



'a i a l -x---a m if j ^ m - + 1 

e>iOi_i ■ • • cr m _i i/ j = to - 1 

0"j-lO"j-2 ' • - 0m if j = i 

X<TiCn+l)(o'i-lO'i) ■ ■ ■ (CmCm + i) i/j = i + 1 



Proof. For x € M let G S„ denote the permutation that x induces on the 
n strands. A so-called permutation braid is a braid in which any two strands 
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cross at most once. If a; is a permutation braid, then x is uniquely determined 
by n x . Moreover, one has Oi ^ x if and only if 7r x (i) > n x (i + 1), that is, if and 
only if strands i and i + 1 cross in x. In particular, the least common multiple 
of two permutation braids x and y can be computed easily from n x and ir y [3] . 
All the braids occurring in the statement of the lemma are permutation braids 
and the claimed equalities are readily checked. □ 

Example 4.8. Consider o 4 o 3 o 2 a 2 ct\ £ B§. From Proposition \4-5\ one has 
F^ im ((X4 l ) = {(7i, 02, 0304}. Using Proposition ^. 6\ and Lemma \4- 7\ repeatedly, 
one then obtains: 

F™ m (o- 4 (7 3 ) = mm {cti, 0203,03^1,03^2,0^0304} = {oi, a 2 a 3 , 04} 

iC m ( fJ 4O-302) = min{ 01 02,02 Vl, 2 \02 O 3 ,02 \o" 4 } = {o-102,03,O- 4 } 

F™ ln (o- 4 a 3 o- 2 a 2 ) = min {0102, 02V102, 02V3, 02V4} = {0102,0302,04} 
F™ ln (04 03 o 2 02 01) = min {oi\0i0 2 , 0i\0302, o"i\o 4 } = {02,0-30201,04} 

Definition 4.9. Let [n] — {1, . . . , n — 1}. 

(a) A function f : [n] —> { — 1, 0}U[n] is called admissible, if it satisfies f(i) < i 
for all i G [n], and /(l) — 1. 

(6) i<br an admissible function f , we define 

F f = {0,0,-1 •••o /w I i e / _1 (M)} u {(Ji-Kii I i e c M. 

We now show that for any w 6 L„, there is an admissible function f w such that 
F™ ln (w) = and we express f( wlT .), for two., e L„, in terms of /„,. 

Proposition 4.10. 

(a) F™ in (e) = F /; where f : [71] ->• {-1,0} U [71] ^zuen &y /(i) = /or al/ i is 
admissible. 

(b) F™ in (0j) = F gj , where g j : [n] -> {-1,0} U [n] given by 

(i ifl<i<j-2 
9j(i) = I -1 ifi = j>l 
I otherwise 

is admissible. 

Proof. Claim (ja|) is trivial. Claim (0) is just a restatement of Proposition 14.51 

□ 
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Proposition 4.11. Let w £ L n and let f be admissible with F™ ln (w) = Ff. 
Then wo~j £ L n if and only if f(j) ^ j. Moreover, in this case, F™ ln (wo~j) = F g , 
for the admissible function g given as follows: 

fori<j-l: g(i)=i 

f/00 iffU)>0 



9(3 - 1) = 
9(3) = 



otherwise 

'0 if j = 1 or f(j) = j 
— 1 otherwise 



iff '(*) = -! 
for i> j: g{i) = { j if f(i) = j + l 

f(i) otherwise 



Proof. One has w<jj £ L n if and only if aj ^ F n (w) ; the latter is the case if and 
only if f(j) ^ j, as o~j is an atom. Now assume wo-j £ L n , and thus /(j) =/= j. 
Then it is clear that g is admissible, as f(j) < j — 1. 

Recall from Proposition 14.61 that F™ ln (waj) is the set of minimal elements in 
{at, . . . , o-j-io-j} U {(TjV/S I /3 £ F™ ln (w)}. This implies that at,..., aj-2 
belong to F™ ln (waj), as they are atoms, so g(i) = i for i < j — 1. Also, 
F™ ln (waj) contains <7j_i<7j, unless j = 1 or crj_i £ F™ ln (waj). We have 
fjj_i £ F™ in (wj) if and only if (Tj\/3 = cr,,_i for some /3 £ F™ in {w). By 
Lemma [4~71 the latter is equivalent to j3 = ajaj-t, and hence to f(j) = j — 1. 
It follows by induction on \w\ that one can have f{i) = —1 only if ai is the 
last letter of w. In that case, at, ■ ■ ■ , Ci_2 forbidden prefixes after w, hence 
j > i — 1. Therefore, if i > J we can have /(i) = —1 only if i = j + 1. As 
(Tj\(TjfTj + i = <7j + i, this is why in this case we have g(i) = i. 
The other values of g follow directly from Proposition 14. 61 using the identities 
from Lemma 14.71 and discarding any elements which are not minimal (specifi- 
cally, any multiples of at, ■ ■ ■ , <Jj-2 and aj-\Uj). □ 

Corollary 4.12. If w £ L n , then there exists a unique admissible function f w , 
such that F™ in (w) = F fa . In particular, \F™ in (w)\ <n-l. 

Proof. If it exists, f w is uniquely determined by F™ ln (w) and thus by w. The 
existence follows by induction on \w\, using Propositions 14. 101 and I4TTT1 □ 

Example 4.13. Consider again a^a^a2a2a\ £ B§. To shorten notation, we 
identify an admissible function f with the sequence [/(l), • • • , f(n — 1)]. From 
Proposition \4-10\ one has / CT4 = [1,2,0,-1]. Repeated application of Proposi- 
tion gTFJ then yields / CT40 . 3 = [1,0,-1,4], U ia3 a 2 = [0,-1,3,4], f ai a 3 ^ 2 a 2 
0. and f ai 

a 3 tT2u 2 cri — [0,2, 1,4]. Observe that the sets of minimal 
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forbidden prefixes described by these functions are exactly those computed in 
Example \4-8\ 




Given w 6 L n , the set F™ ln (w) can be computed efficiently using Proposi- 
tions and the time required is obviously linear in n and \w\. 

There are many admissible functions, but the ones corresponding to minimal 
sets of forbidden prefixes are very special: It can be shown that they are exactly 
those satisfying the conditions in the following corollary. However, proving that 
these conditions are sufficient is quite technical; as we do not use this fact in 
the sequel, we only show that they are necessary. 

Corollary 4.14. Let w = w'o~j £ L n . Let i\ < i-i < ■ ■ ■ < i r be all the indices 
greater than j such that < f w {U) < H (t = 1, • • • ,r). Denote to = f w (j — I )• 

(a) fw(i) = i for i < j - 1. 

(6) f w (i) 6 {0,..-,i} for all i^j. 



(d) If fw{i) = for some i > j, then f w (£) — for £ = i, . . . , n — 1. 

(e) Ifr > 1, then < f w {i r ) <■■■< f w {h) < j. 

if) If r ^ 1 an d m > 0; then either f w (i\) = j or f w {i\) < m. 

(g) Ifr > 2 and m > 0, then < f w (i r ) < ■ ■ • < fwi}2) < m - 

Proof. If w' = e, that is w = aj, the results holds trivially by Proposition 14. 101 
so let w' = w"<Jk and assume the result for f w i. To shorten notation, let 
/ = fw' and g — f w . Since w o~j € L n , we have f(j) j and thus j > k — 1. 
Proposition 14.111 immediately yields Claims (Jaj) , (0) and (Jg) . 

For Claim (0) assume that g(i) = for some i > j. By Proposition 14. Ill this 
happens only if g(i) = f(i) = 0, and if i > k this implies that g(r) = f{r) = 
for r = i, . . . ,n — 1. As i > j > k — 1, the only remaining case is i = k = j + 1 
and f(k) = 0. In this case Claim (jc|) applied to / gives f(k — 1) = k — 1, that 
is f(j) — j, which is a contradiction. Thus Claim (0) holds. 

Now suppose r > 1. We have i\ > j > k — 1, and i\ = k would imply 
f{i\) <E {0,-1} and then, with Proposition 14.111 g(i\) G {0,«i}; the latter is 
a contradiction, so k < i\ < ■ ■ ■ < i r . For t = 1, . . . ,r, Proposition 14.111 also 
yields that either g(i t ) = f(i t ), or f(i t ) = j + 1 and g(i t ) = j. Since f(i t ) < it 
and j < i\ < i^ < ■ ■ ■ < i r , the latter implies that either < f(it) < it, or 
t = 1 and i\ = j + 1 = f(i±) = g{i\) + 1- In any case, applying Claim (je|) to / 
and using k < j + 1, we obtain < f(i r ) < • •• < /(*i) < j + 1- Again using 
Proposition 14. 1 II then yields < g(i r ) < ■ ■ • < g(ii) < j, showing Claim (Je|). 
Suppose to = g(j — 1) > 0. If i\ = j + 1 = /(«i) = + 1, then Claim (0) 
holds. By the preceding paragraph, we may thus assume that < f(it) < U 





if either j = 1 or f w (j - 1) = j - 1 
— 1 otherwise 
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for 1 < t < r. By Proposition 14.111 to > implies f(j) > and to ~ f(j). 
In particular, j ^ k (as f{k) £ {0, —1}), whence we have either j = k — 1 or 
j > k. In the former case, m = f(k — 1), so we can apply the result to / 
and obtain that /(if) < m for 2 < t < r, and either f(ii) — k = j + 1 or 
/(ii) < to; as m < J + 1, Proposition 14.111 then yields g(it) — f(it) < m for 
2 < t < r, and either g(ii) = j or g(i\) = < to, showing Claims (Q) 

and (jg]) in this case. On the other hand, if j > k we just need to notice that 
< m = f(j) < j and apply Claim (pg) to /, which yields f(i t ) < f(J) = m < j, 
whence g(i t ) — f{it) < m. So Claims (0) and (jgj) hold. □ 

One can show that the conditions of Corollary 14.141 are sufficient for / to be 
the defining function of some F™ ln (w) by constructing, for every / satisfying 
the conditions, a word Wf such that F™ ln (wf) = Ff. As this construction for 
general / is very technical, we do not describe it here. Instead, we only consider 
some special cases; enough to show that the number of possible sets F™ ln (w) is 
exponential in n. 

Proposition 4.15. For 1 < i,j < n — 1 define £ A* as the product 

OiO~i±\ ■ ■ ■ o~j . To shorten notation, we identify an admissible function f with 
the sequence [/(l), • • • , f(n — 1)] . 
(a) One has [n — 1, 1] £ L n and 



n—lj] 



[l,...,j-2,0,-l,i + l,...,n-l] ifn-l>j> 1 
[0,2,...,n-l] if j = l 



(b) If i £ [n] and w £ L n such that f w = [0, 2, 3, . . . , i, m.;+i, . . . , m n -i \ , with 
rrik £ {1, k} for k = i + 1, . . . , n — 1, then w ■ [1, i — 1] £ L n and 

[l, . . . , j - 1,0, j, j + 2, . . .,i,m i+1 , . . . ,m„_i] ifl<j<i-l 

[l, . . . ,i - 2,0, i - l,m i+ i, . . . ,to„_i] if j =i-l 



fw[l,j] 



(c) IfiG [n] and w £ L„ such that f w = [l, . . . , i — 2, 0, £ — 1, TOj+i, . . . , to„_i] , 
iwt/i to/j £ {1, fc} /or k = i + 1, . . . , n — 1, then w • [i — 1, 1] £ L n and 

{[!,-•■ , j - 2,0, -1, j + 1, . . . ,i - m i+ i, . . .,m„_i] 
*/ i - 1 > j > 1 
[0, 2, 3, . . . , i - 1, 1, TOj+i, . . . , m„_i] «/ j = 1 

Proof. The claims easily follow by induction on j, using Proposition 14. 1 ll □ 

Corollary 4.16. If S = . . . , i r } is a subset of {1, . . . , n — 2}, possibly empty, 
with i\ > • • • > i r , then ws = [n — 1, 1] • [1, ii] • [ii, 1] • • • • [1, i r ] ■ [i r , 1] £ L n , and 

(0 if j = l 
I j otherwise . 
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Proof. The claim follows from Proposition 14. 1 51 by induction on r. □ 

Corollary 4.17. Let n > 1. The set {F™ in (w) | w G L„} is finite, but it has 
at least 2 n ~ 2 elements. 

Proof. The set is finite as each F^ lin (w) is equal to Ff for some admissible 
function /. It has at least 2™ -2 elements by Corollary 14. 1GI □ 

Before using, in Section [SJ our description of forbidden prefixes to compute the 
numbers x nj k(w,m) required by Algorithm [TJ we show in the next section that 
minimal forbidden prefixes yield a finite state automaton accepting the language 
of lex-representative words, that has the minimal possible number of states. 



5 A minimal finite state automaton accepting L n 

The sets of minimal forbidden prefixes after a given word w G L„ provide a very 
natural way to construct a finite state automaton that accepts the language L n 
of lex-representatives of braids in M. Indeed, we will see that this finite state 
automaton is minimal, in the sense that it has the minimal possible number of 
states. 

Recall that a finite state automaton is a quintuple V = {SS, A, /i, Y, So) where 
SS is a finite set, A is the alphabet, fi : SS x A — > SS is the transition function, 

Y C SS is the set of accepted states and Sq g SS is the initial state [5j. We 
extend /i to a function SS x A* — > SS, also denoted by /i, in the natural way. 

Definition 5.1. Denote T n = (SS,A,^,Y,Sq), where A = {<7i, . . . , cr n _i}, 
SS = {S C A* | S = F™ in {w) for some w e L n } U {A}, S = - i^ in (e) } 

Y = SS\{A}, and \i : SS x A — > SS is given by 

if <Ti € S 

in ({(7i, . . . , o-j_2,o'i_ia-i} U {<r,-\j9 | /? G S 1 }) otherwise . 

Proposition 5.2. T n is a finite state automaton accepting the language L n . 
Moreover, any finite state automaton accepting the language L n has at least as 
many states as T n . 

Proof. By Corollarv l4.12l every state S = F™ ln (w) G SS is uniquely determined 
by an admissible function. As the number of admissible functions is finite, r n 
is indeed a finite state automaton. 

If w G L n , then it follows from Proposition 14.61 by induction on \w\ that 
fj,(So,w) — F™ m (w) G Y, so w is accepted. Conversely, if w ^ L n , let 
w = vtjiv' , where v is the initial subword of maximal length of w that lies 
in L n . As above, we have fj,(So,v) — F™ in (v). Moreover, since vai L n , we 
have Oi G F n (v) and thus at G F™ ln (v) = /j,(Sq,v), as Ui is an atom. Then, by 
the definition of /i, we have fj,(So,vai) — A, and further (by induction on \v'\) 
(j,(Sa,w) = (j,(Sa,v<Jiv') = A, so w is not accepted. 
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In order to show that T n has the minimal possible number of states, assume 
that T' is a, finite state automaton accepting the language L n - Let 1V2 £ L n 
be such that F™ in ( Wl ) ^ F™ in (w 2 ). This implies F n ( Wl ) ^ F n (w 2 ), so by 
symmetry, assume that there exists some a £ F n (wi)\F n (w2) , and consider the 
word w = <jj(a) G L n . As a is forbidden after wi, we have w\w £ L n , so reading 
w starting in the state corresponding to w\ one ends at a fail state. However, 
as a is not forbidden after W2, one has W2W s L n , so reading w starting in 
the state corresponding to W2 one ends at an accepted state; in particular, the 
states corresponding to W\ and u>2 in V must be distinct. As V has at least one 
fail state, the number of states of V is at least the number of states of r„. □ 

Example 5.3. For n = 3, the automaton T n has 5 accepted states; it is rep- 
resented in Figure [H Recall that the initial state corresponds to the empty set, 
and notice that we did not represent the jail state: as usual, arrows which are 
not drawn lead to the fail state {o~\, 02}. 




w 

Figure 2: Finite state automaton accepting lex- representatives in B% . 

Example 5.4. For n = 4, the automaton T n has 18 accepted states; it is 
represented in Figure^ As above, the initial state corresponds to the empty set, 
and we did not represent the fail state {<Ti, 02, 03 }• 

Proposition 5.5. Any finite state automaton accepting the language L n has at 
least 2 n ~ 2 states. 

Proof. By Proposition 15.21 it is sufficient to show that r„ has at least 2™~ 2 
states. The latter follows from Corollary 14.171 □ 

Remark 5.6. We mentioned that the admissible functions that correspond to 
sets of minimal forbidden prefixes are precisely those satisfying the conditions 
of Corollary 14.141 This allows to give an exact expression for the number of 
states of r„. The arguments are extremely technical, however, and as we do 
not use this result in the sequel, we skip the details here. 

Just to give an idea, the following table contains the number of accepted states 
of r n , for n = 3,..., 20. 
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Figure 3: Finite state automaton accepting lex- representatives in Bj. 



n 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


\y\ 


5 


18 


56 


161 


443 


1190 


3156 


8315 


21835 


57246 


149970 



n 


14 


15 


16 


17 


18 


19 


20 


\y\ 


392743 


1028351 


2692416 


7049018 


18454775 


48315461 


126491780 



Remark 5.7. There are generic methods to generate uniformly random words 
of a regular language, such as the recursive method jTQl |B]. These algorithms 
have a precomputation phase, in which a (minimal) acceptor for the language 
in question is computed. 

The results of this section show that the number of states of a minimal acceptor 
for L n is exponential in n. Thus, generating uniformly random words of L n by 
generic language-theoretic methods has a time complexity respectively a space 
complexity that is exponential in n. While efficient for small values of n, such 
approaches are not feasible for larger values of n, as the above figures show. 

6 Counting braids with suitable prefixes 

In this section we finally give a method to compute x n ,k(w, m), where m > 1 if 
w = e, and m > j — 1 if w = w'aj. We assume that we have already computed 
h s ,t for < s < n and < t < k, as well as x n ^ s for < s < k. (See Sectional) 
Recall from Definition 12.21 that x n k(w, m) is the number of lex-representatives 
of length k of the form wv where v G L n does not start with <7j for i < m. 
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As the relations of M are homogeneous, this amounts to counting the braids 
x of length k — \w\ for which wcu(x) G L n and lo(x) does not start with cr, for 
i < m. Knowing the number of braids x of length k — \w\, it suffices to count 
those braids x of length k — \w\ for which wuj(x) £ L n , or u>(x) starts with Ui 
for i < m. If m = 0, the second condition is never satisfied and we simply need 
to count F n (w) n Mfc_| w |. If m > 0, however, we have to add u\, .. . ,a m to the 
set of forbidden prefixes. 

Definition 6.1. Let w £ L n and m G {0, . . . , n — 1}. We define F™ ln (w, m) to 
be the set of minimal elements, with respect to =4, in {a±, . . . , a m } U F™ ln (w). 

Notice that F™ ln (w, m) is obtained from F™ m {w) by removing all elements that 
start with a% , . . . , a m , and adding the elements <j\ , . . . , a m . If m > j — 1, where 
<7j is the last letter of u>, this simplifies the description of the set F™ ln (w,m): 

Proposition 6.2. Let w G L n , let j = 1 if w = e, or w — w'o~j otherwise, and 
let m > j — 1. There is an admissible function f such that F™ ln (w,m) = Ff. 
Moreover, the following hold: 

(a) f(r) = r for r = 1, . . . , m and f(r) > for r = m + 1, . . . , n — 1. 

(b) f(j) =0 ifm = j-l and f(j) =j ifm>j. 

(c) If i £ {j + 1, . . . , n — 1} and f(i) = 0, then f(r) — for r = i, . . . , n — 1. 

(d) J/i G [n] and < f{i) < i, then f(i) < j. 

(e) Ifi,i' G [n], such that < f(i) < i and < /(«') < i', f/ien i' > i implies 

m>m- 

In particular, all elements of F™ ln (w, m) are of the form a r a r ^\ ■ ■ ■ o~fi r ) with 
f(r) < r, and those elements that are longer than one letter are nested: If 
f{i t ) < i t (L <t < s) for j < ii < • ■ • < i s then j > /(i x ) > f(i 2 ) >■■■> f{i s ). 

Proof. This follows immediately from Corollary |4. 141 together with the fact that 
Oi G F™ in (w, m) for i = 1, . . . , m. □ 

We are now in a position to calculate x n ,k{w,m). For a non-empty finite set 
S = {ax, . . . , a r } C M, we define V S = a\ V • • • V a r . We also define V = 1. 
Recall that = {a; G M | |x| = k}; so x„ : fc = \Mk\ and x rii fc = if k < 0. 

Proposition 6.3. Let w G L n and m G {1, . . . , n — 1}. One has 

k—\w\ 

x n>k {w,m)= ^2 (- 1 ) |S ' a; n,fe-|t |-| vs| = ^2 {x n ,k-\w\-i • Ti(w,m)) , 

SC_F» i »(iu,m) 1=0 

where 

T l (w,m)= J2 

SC_F™ in (to,m) 

\ys\=i 
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Proof. The second equality is just a reordering of terms. For a G M, let 
(aM)k = aM n Mfc. As M is a cancellative monoid with homogeneous re- 
lations (aM)k = aM k _\ a \ = {a/3 \ f3 e M fe _| Q |}, that is, |(aM) fe | = x n<k _\ a \. 
The braids x for which wu(x) (fc L n or oj(x) starts with <jj for i < m are exactly 
those that admit a prefix a € F™ ln (w,m). Moreover, = |iu| + |a;|. Hence 



Xn.k (^5 ^"l) — ^n.fc— litil 

By the inclusion-exclusion principle, 



(J (aM) fe _|„| 

a£.F™ 1I1 ('u;,m) 



U M' f )fc-i 

aG-F^ Qin (u',m) 



E 

0^SCF» in (iu,m) 



(-1) 



S -1 



Noting that V = 1 G M, whence | \/ 1 =0, the claim then follows from 



Pi (<xmv h 

a£S 



f] aM)f\M k 



(\/S)M k _ H _ ]ys \ 



((V5)M)nM fc . 



^n,fc-|ra|-| ys| 



□ 



In order to obtain a sub-exponential algorithm, we need to compute Ti(w,m) 
without explicitly summing over all subsets of F^ m (w, m). To this end, we use 
some special properties of the braid monoid. By Proposition ^. 2[ the elements of 
F™ m (w, m) are permutation braids, thus so is the least common multiple of any 
collection of them. The crucial point is that, under certain conditions, we only 
need to know the final position of some strands in a permutation braid, in order 
to describe sufficiently its least common multiple with another permutation 
braid. Specifically, we have the following lemma, which can be readily verified 
working with the induced permutations [5]- 

Lemma 6.4. Let x £ (cr a , o~ a +i 1 ■ ■ • , o~b~i) + C M be a permutation braid and 
let tt x be the permutation induced by x. Let Tr x (a) = a + r and ir x (b) = b — s. 

(a) If a > 1, then x' — a a -i Vi = x(a a -io~ a ■ ■ ■ <T _x+ r ). Moreover, 

b — s if a + r < b — s 

- s 



7iv (a — 1) = a + ■ 



7iv (6) = 



b — s — 1 if a + r > b 
(b) If b < n, then x' — &b V x = x(ab&b-i • ■ ■ Cfc-s)- Moreover, 



7iv (a) = 



if a + r < b — s 



a + r + 1 if a + r > b — s 



ir x '{b + 1) = b — s . 



(c) If b < n then x' = (o~bO~b-i • ■ • o~ a ) Vx = x(o~bO~b—i ■ • ■ o~ a ). Moreover, 
Kx 1 ifl) = a + r + 1 and % x i (6+1) = a . 
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Evaluating Ti(w,m) 

Let w € L n , let j = 1 if w = e, or w = u/<7j otherwise, and let to > j — 1. 
We work with a 3-dimensional array T of size (k — \w\ + l) xnxn, which is trans- 
formed in at most n — 1 steps. At each step, we have two distinguished positions 
a and b with 1 < a < b < n, defining S a , b — F™ m (w, to) D (<r a , cr a +i, • • • , Cb_i) + , 
and for I £ {0, . . . , k — \w\} and r, s <E {0, . . . , n — 1}, the entry T; jI%s of T is 

Tl,r,s = £ (-!) |S| - 

scs a , b , \ ys\=i 

s)(a)=a+r , lr (V S) (&)=&— « 

Initially, we take a = b = j, so = 0. Obviously, in this case all entries of T 
are 0, except T 0j o,o = 1, which corresponds to V = 1. 

At the final step, we will have a = 1 and b = n, so S a .b = F™ un (w, m), whence 

n—1 n—1 

Ti(w, m) = 22 z2 T i,r,s ■ 

In each step we replace (a, b) by (a', 6') € {(a — 1, 6), (a, b + 1)}: If there is some 
o^Ofe-i • • • cr s G F™ ln (w, to) with s < a or if b — n then we set (a', 6') = (a — 1, fe), 
otherwise we set (a', 6') = (a, 6 + 1). By Proposition 16. 21 either S a ' : b'\S a ,b = 0, 
or S a > tb >\S aib = {a b }, or S a > tb >\S at b = {o b o b _ x ■ ■ -cr a }, or S a ' tb >\S atb = {cr a _i}. 
Hence, if S a > tb '\S aib = {x} and ?/ G 5 a .b, then according to Lemma T6.41 the 
length, the displacement of strand a', and the displacement of strand b' of the 
element x V y only depend on x, and the length, the displacement of strand a, 
and the displacement of strand b of the element y. Therefore, the table T can be 
updated when replacing (a, b) by (a! ,b'): The elements contributing to T/ jrs in 
the old table will in the updated table contribute to T/^.s if (a', b') — (a — 1, b), 
respectively to Ti t r,o if (a',b r ) = (a, b + 1). If S a >,b'\S a ,b ^ 0, then in addition, 
their least common multiples with x will contribute to the entry of the updated 
table given by Lemma 16.41 

The arguments from this section yield Algorithm [2] computing x n .k{w,m). 

7 Complexity and timing analysis 

The analysis of the worst-case complexity of Algorithm [T] is relatively straight- 
forward. We assume that the addition of two TV-bit integers has cost O(N), and 
that the multiplication of two Af-bit integers, in the relevant range of N, has 
cost 0(N a ) with a = log 2 3 « 1.585 (Karatsuba multiplication). 

Proposition 7.1. 

(a) The worst-case complexity of Algorithm^ is at most 0(n A k a+1 ) . 

(b) The worst-case complexity of Algorithm\7\is at most 0(n A \nn k a+2 ). 
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Algorithm 2 Computing x n ,k(w,m) 



Input: Integers n > 2 and k > 0. A word w = a ai ■ ■ ■ o at G A* n . An integer m < n—1, 

where m > 1 if w = e and m > j — 1 if tu = wuj. The numbers x„,o, • • • , £?i,fc- 
Output: The number a; ni fc (to, m). 

if £ > k then 
return 

F:=0 

for i := 1 to t do 
if cr ai G F then 

return 
else 

F := min 4 ({ai, . . . , <r ai _ 2 , <7 ai -icx ai } U {cr ai \/3 | /3 € F}) 
F := min^ ({en, . . .,a m } U F) 

Ti, r , s := for all (I, r, s) <E I = {0, . . . , k - t} x {0, . . . , n - 1} x {0, . . . , n - 1} 
if t = then 

a := 1 ; b := 1 ; a := ; T ,o,o := 1 
else 

a :— a t ; b := a t ; a := ; To, 0,0 := 1 
while (a, b) =fc (l,n) do 

if b = n or awh-x ■ ■ ■ cr s G F for some s < a then 
(a', b') := (a - 1, 6) ; a := a + 1 

if F n (cr a /, . . . ,a 6 /_ 1 ) + = Fn (cr a , . . . then 

Eto^Us ifr = 
otherwise 
20: else if (F n (cr a /, . . . , cr 6 /_ 1 ) + ) \ (F n (<r , . . . , <76-i) + ) = {cr a _i} then 

f>- 



1 

2 
3 

4: 

5 
6 
7: 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 



19: for (Z, r, s) G I do T,' r 



21: for (Z,r,s) € I do T/ rs := 



'E:=o^, u , s ifr = 

-T ; _ rjr _i iS if 1 < r < I and r + s < a 

-T;_ rir _i iS _i if 1 < r < I and r + s > a 

, otherwise 



22 

23: 
24: 

25: 



else 

(a', b') := (a, b + 1) ; a := a + 1 

if F n (ov, . . . , ov_i) + = F n (a a , fT6_i)+ then 

for (i,r,8) £ / do r irt := )^=o Tl ^ if S = 

I otherwis 



otherwise 

26: else if (F n {a a ,,. . . , + ) \ (F n (<r , . . . , <r&-i) + ) = {cr 6 } then 

(Y. b -= a Tl,r,u if« = 



27: for (Z, r,s) £ I do T/ 



— Ti-s,r,s-i if 1 < s < Z and r + s < a 
— T ; _ Sjr _i iS _i if 1 < s < Z and r + s > a 
s otherwise 
28: else if (F n (<7 /, . . . , a b '-i) + ) \ (F n (a a , . . . , <76-i> + ) ={a b --- a a } then 

[Etro^.r,™ ifa = 

29: for (Z, r, s) G / do T[^ s := I - r *-a,r-i,« if s = a, Z > a, r > 1 

[ otherwise 

30: (a, 6) :=(o',6'); T := T' 

31: return Ei=o (E"=o E"=o ^V.s) x n,k-t-l 
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Proof. The absolute values of all entries of the array T are bounded by 2™, so 
the entries of T have at most n bits. As it is known (9[ Theorem 10] that the 
logarithmic volume limfc^oo l ° s ^ n - k f B+ is bounded above by 4, the numbers 
x n ,j for < j < k have O(k) bits. 

Consider Algorithm [5] and recall that we can describe forbidden prefixes by 
admissible functions as explained in Sectional By Proposition ^. 11[ line [SI which 
is executed at most k times, and line |H] both have a cost of 0(n). Line ITTJ1 has 
cost 0{n 2 k). Now consider lines [T9l [2TI [251 [27l and [29l At most n executions of 
these occur. Each line involves 0(nk) sums of 0(n) terms, and possibly 0(n 2 k) 
assignments. Since the operands have at most 0(n) bits, the total cost is at 
most 0(n 4 fc). Line |3~T1 involves 0(n 2 k) additions with operands of size 0(n), 
and k multiplications and additions with operands of size 0(n + k), so the cost 
is at most 0(n 3 k + k(n + k) a ) = 0(n 3 k a+1 ). All other lines have cost O(l) and 
are executed at most n times, so Claim (ja|) is shown. 

Consider Algorithm [TJ The body of the for-loop (lines [5] to [T5]) is executed k 
times, the body of the while-loop (lines [71 to [T2l at most O(lnn) times. The cost 
of the for-loop is dominated by the invocations of Algorithm [21 by Claim (Ja|) , 
the cost is at most 0(n 4 Inn k a+2 ). Line[2]has a cost of O(k). From Section[3l 
equation ([3]), calculat ing Yira.i involves ttl additions. We need to calculate hm,j 
for < m < n and < j < k, so overall 0(n 2 k) additions are required; as it 
follows by induction on n that the operands have at most O(nlnn) bits, this 
has a cost of at most 0(n 3 Inn k). Similarly, from equation @, the calcula- 
tion of x n j involves 0(j) additions and multiplications; we need to calculate 
x n j for < j < k, so overall this requires Oik 2 ) additions and multiplica- 
tions of operands with at most 0(nlnn + k) bits, which has a cost of at most 
0((nlnn) Q k a+2 ). This shows Claim ©. □ 

Remark 7.2. There are several optimisations that can be applied when im- 
plementing Algorithms [TJ and [5] While they do not affect the estimate of the 
worst-case complexity, they have a significant impact on actual running times. 

1. It is clearly not necessary to compute the forbidden prefixes from scratch in 
every invocation of Algorithm^ If the forbidden prefixes of the generated 
part of the word are stored in Algorithm [TJ only one invocation of line [S] 
of Algorithm [2] is needed. 

2. It is possible to initialise the array T in Algorithm [2] directly for the values 
of a and b corresponding to the largest range that does not involve a for- 
bidden prefix consisting of more than one letter, using the data computed 
in Section [3] 

3. The array T in Algorithm [5] is in general very sparse, as its contents are 
often collapsed to the 2-dimensional subarrays given by r — respectively 
s = 0. Keeping track of the values of r and s that actually can contain 
non-zero entries and only considering those in any subsequent summations 
greatly reduces actual running times; see Remark 17.31 

4. Since the least common multiple of all forbidden prefixes is a permutation 
braid, its length is bounded above by (")■ Similarly, H n (t) is a polynomial 
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of degree Q), so h n j = if j > Q). It is therefore possible to restrict 
some summations to min{fc, (")} terms. It is clear from the proof of 
Proposition 17.11 that in this case (that is, n constant and k > (")) the 
complexity in terms of k is reduced to 0(k a+1 ). 

Table [1] contains timing results for several values of n and k, using the imple- 
mentation of the algorithm in the C-kernel of the computational algebra package 
Magma [3] by the first author. Computations were done with a development 
version of Magma V2.16 on a GNU /Linux system with an Intel E8400 64-bit 
CPU (core: 3 GHz, FSB: 1333 MHz) and a main memory bandwidth of 4.7 GB/s 
(X38 chipset, dual channel DDR2 RAM, memory bus: 1066MHz). The mini- 
mum sample size was 100. The sample size was varied to achieve an approximate 
running time of 2.5 minutes per sample (where possible honouring the minimum 
sample size), and hence a sufficient degree of accuracy for all parameter values. 





4 


8 


16 


n 

32 


64 


128 


256 


512 


4 


1.09e0 


1.33e0 


1.76e0 


2.07e0 


2.45e0 


2.95e0 


3.58e0 


8.03e0 


8 


4.62e0 


6.39e0 


5.79e0 


6.10e0 


1.26el 


2.19el 


3.13el 


4.84el 


16 


2.20el 


5.10el 


4.56el 


7.52el 


9.73el 


1.35e2 


2.31e2 


2.42e2 


32 


1.07e2 


3.81e2 


4.63e2 


6.70e2 


8.81e2 


1.01e3 


1.15e3 


1.79e3 


k 64 


4.84e2 


2.07e3 


4.55e3 


6.72e3 


5.72e3 


7.94e3 


7.61e3 


7.72e3 


128 


2.01e3 


9.55e3 


3.30e4 


5.24e4 


5.70e4 


1.03e5 


1.13e5 


1.04e5 


256 


8.42e3 


4.27e4 


1.90e5 


4.17e5 


6.06e5 


1.20e6 


1.81e6 


2.06e6 


512 


4.36e4 


1.78e5 


9.70e5 


3.27e6 


5.85e6 


1.23e7 


2.56e7 


3.29e7 


1024 


2.52e5 


7.24e5 


4.77e6 


2.07e7 


5.22e7 


1.25e8 


2.78e8 


4.29e8 



Table 1: Running times in /j,s per random braid for different values of n and k. 



k 


4 


8 16 


32 


64 


128 


256 


512 


1024 




0.35 


0.50 0.48 


0.48 


0.47 


0.75 


1.08 


1.36 


1.59 



Table 2: Observed exponent e„ of n for different values of k; see Remark I7T31 



n 


4 8 


16 


32 


64 


128 


256 


512 




2.30 2.41 


2.80 


3.02 


3.08 


3.19 


3.27 


3.21 



Table 3: Observed exponent e& of k for different values of n; see Remark 17.31 

Remark 7.3. Assuming that the observed running time t is approximated by 
a relation of the form t — c n e ™ k ek , the value of e n for a fixed value of k can be 
obtained by applying regression analysis to the rows of Table [T] Similarly, the 
value of ek for a fixed value of n can be obtained from the columns of Table [TJ 
The observed values of e„ and et are given in Tables [2] and [31 respectively. 
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1. The observed values of ej, are close to a + 2 w 3.585 for large values of n. 
For small values of n, where k > (™) f° r most of the parameter range 
tested, the observed values of are closer to a + 1 » 2.585. 

This is exactly what should be expected from Proposition 17.11 and Re- 
mark [7T2] (|4]). That is, as far as k is concerned, the average case complexity 
is quite close to the worst-case complexity. 

2. The observed values of e n are significantly smaller than 4, that is, the 
average case complexity in n is significantly better than the worst-case 
bound of Proposition 17.11 

This is not surprising in the light of Remark 17.21 ([3]) : Our implementation 
exploits the fact that the array T is sparse and tends to contain non-zero 
entries only for relatively few values of the indices r and s. The experimen- 
tal results suggest that the array T effectively is far from 3-dimensional, 
especially if k is relatively small. We expect that the theoretical growth 
rate in n will be achieved, if at all, only for extremely large values of k. 
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